Bogoliubov shadow bands in the normal state 
of superconducting systems with strong pair fluctuations 
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On the basis of a scenario where electron pairing is induced by resonant two-particle scatter- 
ing (the Boson Fermion model), we show how precursors of the superconducting state - in form 
of overdamped Bogoliubov modes - emerge in the normal state upon approaching the transition 
temperature from above. This result is obtained by a renormalization technique based on continu- 
ous unitary transformations (the flow equations), projecting out the coherent contributions in the 
electron spectral function from an incoherent background. 
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The discovery of high temperature superconductivity 
(HTS) and the enormous theoretical effort following this 
discovery has led to reconsider the whole issue of su- 
perconductivity and to put it into a more general and 
broader perspective than that of the familiar BCS sce- 
nario. It became clear early on that in these new super- 
conductors one is confronted with a situation between 
that of classical Cooper pairing (controlled by the am- 
plitude of the order parameter) and that of a superfiuid 
phase of tightly bound electron pairs (controlled by the 
phase fluctuation). The theoretical issue is to under- 
stand the interplay between the two, which amounts to 
fully taking into account the interaction between single 
electron states and electron pair fluctuations. A partic- 
ularly clear presentation of this problematic was given 
by Tchernyshyov Q], who introduced an effective ampli- 
tude for electron pairs in the normal state by includ- 
ing in the propagator of an electron a process involving 
the simultaneous propagation of a hole together with a 
Cooperon. This leads to the emergence of two-electron 
resonant states (long lived pair fluctuations) inside the 
Fermi sea, manifest in form of a pseudogap - a precursor 
of the superconducting gap in the normal state. 

The phenomenon of superconducting correlations 
building up above T c - the onset temperature of super- 
conductivity - is by no means restricted to the HTS. It is 
encountered in a great variety of physical systems, where 
resonance pairing of fermionic quasi particle plays a key 
role and which can be discussed on the basis of the phe- 
nomenological Boson Fermion model (BFM) . This model 
describes the formation of resonance bound Fermion- 
pairs inside the Fermi sea of uncorrelated Fermions. The 
resonant bound states are monitored by an exchange in- 
teraction with localized Bosons (tightly bound Fermion 
pairs, the specific origin of which will be of no importance 
for the present study). Realizations of such physics are 
be found in: i) electron-phonon coupled systems in the 
intermediary coupling limit (see for instance ref.|2| for 
which this BFM was devised initially) describing an ex- 
change interaction between localized bipolarons and itin- 



erant electrons, ii) the positive U Hubbard modelQ de- 
scribing an exchange interaction between spinon singlets 
of RVB electron pairs and holons, iii) low density nuclear 
matter with isospin singlet pairing i4 and iv) the Fesh- 
bach resonance in atomic physics |5|, [g] invoked for tuning 
a quasi-bound state of atoms through a threshold result- 
ing in a resonance superfluidity in traps 0,01 > and which 
involves entanglements of atoms in squeezed states^. 

A salient feature of this BFM is the appearance of a 
pseudogap in the normal state single Fermion density 
of states when decreasing the temperature below a cer- 
tain value T* . Accompanied with this is a qualitative 
change in the transport: going from single Fermion trans- 
port above T* (describing a relatively good metal) to 
one involving well defined itinerant bosonic charge car- 
riers (describin g a bosonic metal coexisting with a poor 
fermionic metal[lfj,|ll|). Below a certain critical temper- 
ature T c , the itinerant Bosons together with finite mo- 
mentum Cooperons, induced in the fermionic subsystem, 
end up in a phase correlated superfiuid state. One hence 
should not only expect a pseudogap in the Fermionic den- 
sity of states as a precursor to the superconducting phase, 
but also corresponding modifications of the sizable inco- 
herent component of the Fermionic spectral functions, 
which characterizes the pseudogap phase in the temper- 
ature interval [T C ,T*], to evolve into damped, but never- 
theless clearly discernible, Bogoliubov modes as one ap- 
proaches T c . The feasibility of such a scenario will be 
the issue of the present study. 

Applying a flow equation renormalization technique to 
the BFM permits us to construct renormalized Fermionic 
operators which posses features of potentially containing 
residues of Bogoliubov type modes in the normal state. 
The essence of this technique consists in devising an infi- 
nite series of continuous unitary transformations |l 2l fl3j | 
for the Hamiltonian, such as H(l) = e s (^He~ s "\ where 
I denotes the continuous flow parameter. Imposing a 
constraint structure on the renormalized Hamiltonian de- 
fines renormalization equations for the parameters of this 
Hamiltonian which are devised in such a way that the cx- 
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change interaction between the Bosons and the Fermions 
is renormalized to zero giving rise to decoupled systems of 
Bosons and itinerant Fermions, whose parameters how- 
ever have been rendered interdependent in the course 
of such a renormalization procedure. In the momen- 
tum representation this Hamiltonian consists of two parts 
H(l) = Ho(l) + Hi n t (I) and has the following structure 
in the l-th step of this continuous transformation 



H (l) = 5>k(0 - M)cL<** + 52(E«(l) - 

k.cr q 



Hintyl) 



c k2 re f er to annihilation (creation) operators for the itin- 
erant Fermions with the energy £k(0 and 6 q (& q ) denote 
boson operators representing bound Fermion pairs with 
energy -Eq(Z) — 2/i. The Boson- Fermion exchange cou- 
pling is denoted by Wk,p(0 and the interaction between 
Fermions by t/k, P . q (0- For I = these parameters re- 
duce to the bare quantities £k — fi, As — 2fi, v and 
respectively, which characterize the initial Hamiltonian. 

In some recent work^4| we investigated the struc- 
ture of the final renormalized Hamiltonian. Its result- 
ing energy spectrum and single Fermion density of states 
showed the opening up of a pseudogap below a certain 
T* of the order of v. We here shall investigate the spec- 
tral function of these Fermions and how it changes as 
we go from the pseudogap phase [T c , T*] into the super- 
conducting one below T c . We propose for that purpose 
a procedure for evaluating correlation functions within 
such a flow equation technique, which closely follows the 
standard procedure for renormalizing the Hamiltonian, 
determined by the differential equation: 



k,p.q 

-7= Wk ,p(0 ( 6 P+k c kiC P T + k' 
k,p 



(2) 



dH(l)/dl = [?/(/), H(l)], 



(3) 



and subject to the initial condition H(0), presenting the 
original Hamiltonian. A suitably generating operator 
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where a k , P (0 = kk(0 + £ P (0 - #k+p(0) u k, P (0 is cho- 
sen in such a wav[l2j that lim^oo Hi nt (l) = 0. This leads 
to a set of differential equations (given by eqs. (16-21) 
in Ref. which determines the evolution of the I de- 
pendent parameters of the Hamiltonian. For consistency 
reasons, a flow equation procedure for an arbitrary oper- 
ator 0(1) ought to be controlled by a formally equivalent 
equation of that defining the evolution of the Hamilto- 
nian, i.e., eq. El where H(l) is replaced by 0(1). This 
leads to a flow parameter dependent parameterization of 
the various operators for which we impose (similar to the 



procedure for deriving a renormalized Hamiltonian) the 
following constraint structure for the fermion operators: 

Ckt(0 
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which generalizes the standard Bogoliubov transforma- 
tion in two ways: (a) the initial particle and hole opera- 
tors are transformed into Bogoliubov modes in a contin- 
uous way, (b) the correlated motion involving Fermion 
holes and Cooperons is taken into account via the terms 
proportional to Pk,q(0 and r^ q(l). 

These flow dependent parameters are then determined 
by the following set of differential equations: 

dPk(l) 



dl 
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dK k (l) 

dl 
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dpk,q(l) 

dl 

rfr k ,q(0 
dl 
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(8) 
(9) 

with the initial conditions 'Pk(O) = 1, 7£k(0) = 0, 
Pk,q(0) = 0, rk, q (0) = 0. n^ ond denotes the fraction of the 
condensed bosons and ro q = (6 q 6 q ) = [exp(i? q /fcsr) — 
1] _1 the distribution of the finite momentum (q 7^ 0) 
bosons. This set of equations has to be solved in con- 
junction with those (eqs 16-21 in ref. 14) determining the 
evolution of the various parameters entering the Hamil- 
tonian (£k(0i H q (l) and £7k,p,q(0) being linked together 
via the expression for ctk. p (l). The above flow equations 
iJBEJ satisfy the sum rule 
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relating the weights of the various coherent and incoher- 
ent contributions to the spectral function. With these 
definitions, the single Fermion spectral function becomes 

A F (\t^) = |Pk(oo)| 2 «5( W -£k) + |^k(oo)| 2 5( W + e_ k ) 

]V H ( n q + n q+kt) bk, q (oo)| 2 (5(w + £ q - e q+k ) 
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|r k , q M| 2 ^-S q + eVk) (11) 
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which is composed of a coherent part A F oh {k.,w), repre- 
sented by a (5-function like peak, and a remaining incoher- 
ent background A[ nc (k), given by the last two terms in 
cq. (|llfl . Ek and Eq refer to the end results of the renor- 
malization procedure (I = oo) for the Fermion and, re- 
spectively, Boson spectra. We now analyze the structure 
of i|ll|) for different characteristic temperature regimes: 
T > T*, the pseudogap regime [T C ,T*] and the super- 
conducting regime T < T c . 

In the normal phase we have n^ ond = 0, which implies 
7^k(0 = Pk,q(0 = for any I and hence 

Af o/l (k,c) = |7> k M| 2 <5(^- £ ~ k ) (12) 

A( nc (K^) = ^EK+<-ki)w»)i 2 

x 5{u - £ q + £ q _ k ) (13) 

for the coherent and incoherent contributions to the 
Fermion spectral function. As discussed in ref. ^J, f° r 
temperatures above a certain T* the Fermion dispersion 
is essentially unrenormalized and given by exclusively 
A F oh (k,uj) with — 1. Below T* , on the contrary, the 
Fermion renormalization is becoming increasingly impor- 
tant and gives rise to the opening of a pseudogap in the 
single particle density of states with a spectral weight 
becoming smaller than unity and a consequent redistri- 
bution of this lacking spectral weight into an incoherent 
part given by A[ nc (k, uj). The fully selfconsistent solu- 
tions to this problem by numerical means are presented 
in Fig^ Throughout this work we take as energy unit 
the Fermion bandwidth 2zt, {z denoting the coordina- 
tion number) and the initial parameters Ab = —0.6 and 
v = 0.1. We, moreover, choose the total particle con- 
centration per site n tot = n F + 2n B = 1.0. Since we 
are interested here only in the qualitative aspects of the 
underlying physics and since the physics in the normal 
state is controlled by very local correlations and hence 
little dependent on the dimensionality, we can approxi- 
mate the various integrals in this temperature regime by 
their sum over a ID Brillouin zone with 2000 wavevec- 
tors, k n — n^j^ and —1000 < n < 1000 as long as we 
don't get too close to T c . 

In the superconducting phase, from inspection of eq. 
pi|l. we notice the appearance of two contributions to 
the coherent part of the spectral function with spec- 
tral weights 'Pk(oo) 2 and TZ^oo) 2 , with a total weight 
'Pkioo) 2 + 7?-k(oo) 2 amounting to less than unity be- 
cause of the sum rule, Eq. I|10l) . Well below T c the 
quasi particle energies are given by[14| £k = sgn(ek — 

f^)\J ( £ k — £t) 2 + v2n cond an d present the standard Bo- 
goliubov modes which follow from a straight forward 
mean field analysis of this BFM. The characteristic fea- 
tures of the Fermion spectral function in the supercon- 
ducting phase can be assessed form their zero temper- 
ature limit. There, n B — for any q / and the 
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FIG. 1: The single particle fermion spectral function 
A F (k, uj) decomposed into its coherent (thick bars whose 
height indicate the intensity of the delta like contributions) 
sitting on top of an incoherent component in the vicinity of 
kF (indicated by the bold spectral line) for the normal phase 
(a) above T* (T = 0.02), (b) and (c)for the pseudogap region 
T* > T > T c (0.007, 0.004) and (d) for the superconducting 
phase (in the ground state T = 0). The distance between 
the neighboring lines corresponds to changes in wavevector 
by multiples of Ak = 7r/1000a. 
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Fermion distribution function reduces to a step function 
n£ a = 0(-£k)- Hence, 

< c (k,u,) = l£ [b k , q M| 2 ^ + S q -£ q _ k M-£k- q ) 

+ kk,q(00)| 2 5(0> - -Bq + £q+k)6'(-£q+k) (14) 

In the spirit of the qualitative study, presented here, the 
sum over the wave vectors q ^ in this expression for 
Af nc (k,uj) can again be safely approximated by a sum 
over a, ID Brillouin zone. 

In Figs, la- Id we illustrate the evolution with tem- 
perature T of the coherent (the (5-function like peaks) 
and the incoherent components (the broad hump like 
features) of the Fermion spectral function in the rele- 
vant wave vector regime around kp where qualitative 
changes are manifest. In the high temperature regime 
(T > T* , Fig. la) the single Fermion spectrum is practi- 
cally unrenormalized, but is constraint to the minimum 
of a broad "high temperature incoherent contribution" 
which is visible for energies above as well as below this 
coherent contribution. As we go, upon lowering T, into 
the pseudogap phase (T c < T < T* , Figs, lb, and lc), 
an additional "low temperature incoherent component" 
emerges, having a dispersion opposite to that of the "high 
temperature incoherent contribution" . That latter re- 
mains practically unchanged all the way down in tem- 
perature, and right into the superconducting phase. The 
"low temperature incoherent component", on the con- 
trary, noticeably narrows upon decreasing T, and eventu- 
ally strongly modifies the behavior of the coherent com- 
ponent, foreshadowing Bogoliubov modes of the super- 
conducting phase. The spectral shape of this low temper- 
ature incoherent component" becomes increasingly bet- 
ter defined, as we approach the superconducting phase. 
Below the pseudogap it shows a dispersion which bends 
downwards for increasing wavectors above kp. Above 
the pseudogap this dispersion of the shadow Bogoliubov 
mode bends upwards for decreasing wavectors below kp. 
These features are clearly apparent from a comparison of 
Fig. lc with Fig. Id for the superconducting phase at 
T = 0. In Fig. Id the standard well defined Bogoliubov 
modes are clearly visible. 

In order to highlight this emergence of the Bogoliubov 
modes out of the incoherent contribution of the Fermion 
spectral function in the normal state we illustrate in 
Fig. 2 the Fermion spectral function at kp for different 
temperatures. We notice the gradual narrowing of this 
low temperature incoherent contribution to this spectral 
function as we approach the superconducting phase upon 
lowering the temperature, with peak positions practically 
being independent on temperature. 

In conclusion, we have shown that in systems with pre- jg 
cursor pairing, such as given by the Boson Fermion sce- 
nario, we can expect rather well defined remnants of the [10 




FIG. 2: Evolution with temperature of the Fermion spectral 
function A F (k^ , co) at kF in the vicinity of the Fermi energy. 
The spectral weights of the 5-function like peaks (correspond- 
ing to the coherent components sitting on top of the incoher- 
ent ones) are indicated by squares, circles and triangles. 



Bogoliubov modes in the normal phase in a restricted 
regime above T c . These modes eventually broaden into 
intrinsically incoherent contributions of the spectral func- 
tion as the temperature increases and approaches T* 
where resonance pairing of the Fermions ceases. Experi- 
mental verification of such shadow Bogoliubov bands in 
the high HTS would be decisive in determining whether 
the pseudogap and the superconducting gap in these ma- 
terials are of common nature or not. In the present calcu- 
lation we neglected the effect coming from the two-body 
interaction C/k.pjjXoo) whose magnitude is small and of 
the order ~ v 2 14]. Nevertheless in the pseudogap fre- 
quency regime this interaction would lead to a shifting 
away of the quasiparticle peaks from this frequency re- 
gion (and hence reenforce the pseudogap features quali- 
tatively) and simultaneously lead to a broadening of the 
delta peak structure of the coherent contributions of or- 
der of the Cooperons bandwidth ~ v 2 . This will be dis- 
cussed in detail in some future study. 
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